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Abstract. In this paper we discuss the X-ray properties of 49 local {z < 0.035) Seyfert 2 galaxies with HST/WFC2 
high-resolution optical coverage. It includes the results of 26 still unpublished Chandra and XMM-Newton obser- 
vations, which yield 25 (22) new X-ray detections in the 0.5-2 keV (2-10 keV) energy band. Our sample covers a 
range m the 2-10 keV observed flux from 3 x 10~" to 6 x 10"^^ erg cm s . The percentage of the objects which 
are likely obscured by Compton-thick matter (column density, Nh > Cj"^ ~ 1.6 x 10^^ cm~^) is ~50%, and reaches 
~80% for log(F2_io) < 12.3. Hence, fluorescent iron lines with large Equivalent Width {EW > 0.6 keV) are 
common in our sample (6 new detections at a confidence level >2a). They are explained as due to reflection off 
the illuminated side of optically thick material. We confirm a correlation between the presence of a ~100-pc scale 
nuclear dust in the WFC2 images and Compton-thin obscuration. We interpret this correlation as due to the 
large covering fraction of gas associated with the dust lanes (following an idea originally proposed by Malkan et 
al. 1998, and Matt 2000). The X-ray spectra of highly obscured AGN invariably present a prominent soft excess 
emission above the extrapolation of the hard X-ray component. This soft component can account for a very large 
fraction of the overall X-ray energy budget. As this component is generally unobscured - and therefore likely 
produced in extended gas structures - it may lead to a severe underestimation of the nuclear obscuration in 2 ~ 1 
absorbed AGN, if standard X-ray colors are used to classify them. As a by-product of our study, we report the 
discovery of a soft X-ray, luminous (~7xl0'"' erg s~^) halo embedding the interacting galaxy pair Mkn 266. 
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^ 1. Introduction on that scale. However, both the low covering fraction 

! of the BLR and their optical thickness seem too low to 

; While it is widely accepted that "type 2" Active Galactic account fo r the large fraction of heavily X-ray obscured 

^ . Nuclei (AGN) are generally obscured in X-rays by large aGN (e.g. IRisaliti et al. TMl . Other pieces of evidence 

■ - - ' column densities, Nh, the exact location and distribu- _ ^gain mainly based on X-ray variability studies - sug- 

tion of the obscuring gas are still debated. Traditional ggg^ ^ complex structure for the X-ray obscuring matter 

wisdom identifies the obscuring gas with matter responsi- j^any instances, and in particular rules out a single 

ble for orientation-dependent effects {e.g., the dusty com- homogeneous, azimuthally-symmetric obscuring structure 

pact molecular "torus" of the AGN Unification Scenarios; (|Matt et al. 2008al K^uainazzi et al. 2005all. A lthough an 

[Antonucci fc Miller 1985. Antonucci 1993). However, re- inhomogeneous "patchy" torus llPier fc Krohk 1992.1 is 

cent observational results challenge this interpretation. ^^^^ consistent with the data, and with the pres- 

Fast changes in the X-ray column density advocate rather ence of Compton-t hick reflection in unobscured AGN 

for matter closer to the nuclear superm assive black hole iHjjanclii et al. 200411 . a possible scenario is that X-ray 

such as the Broad Line Region (BLR; [ Elvis et al. 20 01 obscuration occurs on various scales, sometimes due to 

[Lamer et al. 2003t- "Narrow" components of the K„ flu- j^ore compact matter close to the black hole. 



some- 



orescent iron lines in unobscured AGN - when resolved times due to matter associated with the host galaxy 

by high-resolution spectroscopic Chandra measurements ^.^ther then with the nuclear environment. In a few cases, 

- are so metimes consistent with an origin in the BLR "multiple absorb ers" are directly measured in X-rays 

as wen ( [Yaqoob et al. 2001[ 2004; Bianchi et al. 2003b), ([Matt et al. 200ll 2003a), and their lowest column densi- 

providing evidence for the existence of obscuring matter tjgg ^re consistent with the optical reddening due to host 
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galaxy dust lanes, once a standard Galactic gas-to-dust 
ratio is assumed^. 

Another possible approach to elucidate the location 
and geometrical distribution of X-ray obscuring matter 
in AGN is a comparison between the distribution of X- 
ray column densities and indicators of reddening and/or 
extinction at other wavelengths. Prior to the launch of 
the major X-ray observatories Chandra and XMM-Newton 
this task was hampered by the limited number of objects, 
for which X-ray spectroscopic measurements were avail- 
able. 

Guainazzi et al. (2001; GOl hereafter) studied the 
correlation between X-ray obscuration and the presence 
of nuclear dust in the innermost ~100 pc around the 
nucleus, as detected in the high-resolution WFC2/HST 
imaging survey of nearby Seyfert galaxies after Malkan et 
al. (1998; M98 hereafter). GOl discovered a possible cor- 
relation: Compton-thin objects tend to be preferentially 
located in galaxies with a larger dust content, whereas the 
likelihood that a galaxy hosts a Compton-thick nucleus is 
not affected by its nuclear dust content. However, the sta- 
tistical significance of the correlation in GOl was at the 
l-fT level. Since then, we have carried out a XMM-Newton 
observing program to verify this tentative correlation on 
a larger sample. The results of this program, coupled with 
the analysis of still unpublished archival data of objects 
belonging to the M98 sample, are the main subject of this 
paper. 

The structure of the paper is as follows: in Sect. 2 
we introduce the sample; in Sect. 3 we present the reduc- 
tion/analysis procedure. The spectral results are described 
in Sect. 4, and discussed in Sect. 5. Our main conclusions 
are summarized in Sect. 6. 

In this paper: the acceptance level for a spectral fit cor- 
responds to a null hypothesis probability < 5%; errors on 
the spectral fit parameters and thereof derived quantities 
are at the 90% confidence level for 1 interesting param- 
eter; all other uncertainties are at the l-cr level; energies 
refer to the source rest frame, unless otherwise specified. 

2. The sample 

The sample is made of all the Seyfert 2 galaxies in the 
Veron catalogue jVeron-Getty fc Veron 2003| ), which: a) 
have z < 0.035; b) have been observed in X-rays; c) 
have been included in the high-resolution HST/WFC2 sur- 
vey of M98. The plate scale of WFC2 images (0".046) 
corresponds to a pixel size ^ 32 pc, and the total im- 
age size (37"x37") is ^ 25 kpc. The total number of 
objects fulfilling all three criteria is 49^, which repre- 

^ On a possible systematic difference in the gas-to-dust 
ratio between our Galaxy and luminous nearby AGN, see 
iGranato et al. 1997i and references therein 

^ We have excluded from the sample UGC 4203 (the 
"Phoenix Galaxy"), as it belongs to the class of "Changing- 
look" Seyfert 2 galaxies (IMatt et al. 2003bl . and its classifica- 
tion with respect to the X-ray obscuration is therefore ambigu- 
ous 



sent ~16% of the total parent sample. 20 of them were 
discussed by GOl. Although updated measurements for 
the GOl sample have been used, whenever published 
in the literature, we did not systematically reanalyze 
all GOl sources, as they are mostly well studied, X- 
ray bright AGN. Their X-ray absorption classification - 
originally derived from Bassani et al. (1999) - is there- 
fore robust. The X-ray spectra of three more M98 ob- 
jects have been published elsewhere: ESO 509-IG066 
l|(;uainazzi et al. 2005hll. NGC 4 24 IMatt et aJ. 2003^11. 
NGC 7130 ULevenson et al. 2005|l . We present in Sect. 4 
the results of C/ianrfra/XMM-Newton observations of the 
26 still unpublished objects belonging to the sample. 

Following the classification in M98 we consider "nu- 
clear dusty" those galaxies, whose nuclear dust morphol- 
ogy falls into the "DC" ("Dust passing close or cross- 
ing the Center") and "DI" ("Dust Irregular")^ categories. 
We have verified that the adopted classification is consis- 
tent with that reported by other authors, which discuss 
the same nuclear images, or images with comparable spa- 
tial resolution l|Ferruit et al. 20001 IQuillen et al. 1999b[ 
IMartini et al. 2003|l . The mean redshifts for "dusty" and 
"non-dusty" objects are consistent: (z) = 0.016 ±0.009 for 
the former, (z) = 0.017 ±0.008 for the latter, respectively. 
The overwhelming majority of the galaxies in our sam- 
ple (86%) are spirals. The remaining are almost equally 
shared between ellipticals (5%) and irregulars (9%). 

The list of objects, for which new X-ray data are pre- 
sented in this paper, is reported in Tab.^ with redshifts. 
Ha/ Hp flux ratios, 0[iii] fluxes and some information on 
the C/iandra/XMM-Newton observation. 

3. Observations and data reduction 

The XMM-Newton data were reduced with SAS v6.1 
l|Gabriel et al. 2003jl . using the most updated calibration 
files available at the moment the data reduction was per- 
formed. In this paper, only data from the EPIC cam- 
eras (MOS; Turner et al. 2001; pn, Striider et al. 2001) 
will be discussed. Event lists from the two MOS cameras 
were merged before accumulation of any scientific prod- 
ucts. Single to double (quadruple) events were used to 
accumulate pn (MOS) spectra. High-background particle 
flares were removed, by applying fixed thresholds on the 
single-event, _E > 10 kcV, At 10 s light curves. These 
thresholds, as well as the radius of the source circular ex- 
traction regions, were optimized to maximize the signal- 
to-noise ratio, and are listed in Tab. ^ The background 
was extracted from annuli around the source for the MOS, 
and from circular regions in the same chip for the pn, at 
the same height in detector coordinates as the source lo- 
cation. Spectra were binned in order to oversample the 
intrinsic instrumental energy resolution by a factor ~3. If 
the resulting spectral bins had less than 25 background- 
subtracted counts, the spectra were further rebinned to 
reach this limit, which ensures the applicability of the 

^ i.e., "not associated with any spiral arm pattern", M98 
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Table 1. Properties of the Seyfert 2 galaxies discussed in this paper. The "Date" indicates the observation date. Tesp 
is the net exposure time in the pn (ACIS) XMM- Newton (Chandra) camera. The "Threshold" was applied on the 
XMM-Newton EPIC E> 10 keV, field-of-view, single-event light curves employed to remove intervals of high particle 
background. "Radius" refers to the size of the source extraction region in MOS/pn (XMM-Newton observations) and 
ACIS-S (C/iandra observations), respectively. 



Source 




z 


Nn.Gai 




Vo[iii\°' 


Ref." 


Date 


Texp 


Threshold 


Radius 








(10 cm ) 








(DD-MM-YY) 


(ks) 


(s ) 


(arcsecs) 


beyrert / exnibiting 


nuclear dust 


















Mkn266 




0.0281 


1.7 


5.89 


0.32 


8 


02-11-01 


19.7 




2.0 


Mkn612 




0.0203 


5.0 


6.61 


1.80 


8 


19-01-04 


9.1 


0.35/1.0 


40/25 


NGC4156 




0.0225 


2.0 


33.33 




5 


22-12-00 


51.2 


.../... 


40/40 


NGC4968 




0.0099 


9.2 


14.9 


1.00 


8 


05-01-01 


4.3 


.../... 


40/40 


NGC7212 




0.0266 


5.5 


5.01 


7.08 


6 


20-05-04 


10.4 


.../l.O 


30/60 


NGC7465 




0.0066 


6.0 


5.34 




1 


c\r\ /^o 

29-07-03 


5.0 




6.0 


UGC2456 




0.0120 


12.1 


8.51 


2.40 


8 


i^/^ j^i^ r\f^ 

20-02-05 


11.7 


0.35/2.0 


35/35 


UGC4229 




0.0232 


5.2 


18.20 


0.40 


6 


05-05-03 


4.8 


0.35/1.0 


20/40 


UGC6527 




0.0274 


1.1 


6.57 


5.40 


6 


04-11-04 


7.8 


0.35/1.0 


50/40 


UGC987 




0.0155 


5.7 


5.38 


0.30 


2 


23-01-04 


14.6 


0.35/1.0 


50/40 


Seyfert 2 not exhibiting nuclear dust 


















ESO269-G012 




0.0165 


13.8 




0.50 


7 


06-01-04 


17.6 




6.8 


2MASSJ01500266-0725482 


0.0177 


2.3 


7.62 


0.53 


4 


21-01-04 


9.0 


0.35/1.0 


45/45 


IC 3639 




0.0109 


5.1 


6.14 


3.14 


4 


07-03-04 


8.7 




5.2 


IC4995 




0.0161 


4.2 


5.37 


2.81 


8 


25-09-04 


6.7 


0.35/1.0 


20/30 


IRAS 15480-0344 




0.0303 


9.4 


10.19 


1.38 


4 


31-07-03 


2.5 


0.35/2.0 


20/25 


NGC1410 




0.0253 






u.ou 


s 
o 




7 R 


\j.oo 1 ^.yj 




Mknl 




0.0154 


5.5 


3.32 


3.47 


4 


09-01-04 


9.1 


.../2.0 


45/30 


NGC17 




0.0198 


2.6 


43.75 


0.11 


8 


22-12-02 


10.6 


0.35/1.0 


40/40 


UGC1214 




0.0172 


3.0 


4.17 


15.85 


8 


15-01-04 


9.1 


.../... 


32/32 


UGC2608 




0.0233 


14.1 


6.31 


2.19 


6 


28-01-02 


1.1 


1.0/2.0 


25/25 


NGC5283 




0.0104 


1.8 


3.80 


2.69 


8 


24-11-03 


8.9 




2.9 


NGC5427 




0.0087 


2.2 








26-03-04 


8.8 




C 


NGC5953 




0.0066 


3.3 


22.0 


8.60 


5,9 


29-12-02 


9.9 




2.5 


UM625 




0.0250 


4.0 


9.17 


2.10 


3 


02-07-04 


7.4 


.../l.O 


30/30 


NGC3982 




0.0037 


1.2 


4.00 


1.51 


8 


03-01-04 


9.2 




3.7 


NGC591 




0.0152 


4.8 


6.02 


2.30 


6 


12-01-04 


8.8 


0.35/1.0 


25/40 


"in units of 10"^^ ere 


' cm ^ s 


-1 



















'l.- Osterbrock et al. (1987); 2.- Dahari et al. (1988); 3.- Terlevich et al. (1991); 4.- de Grijp et al. (1992); 5.- Kennicutt et al. 
(1992); 6.- Mulchaey et al. (1994); 7.- Mulchaey et al. (1996); 8.- PoUetta et al. (1996); 9.- Risaliti et al. (1999) 
'^not detected 



statistics to evaluate the quality of the spectral fitting, pn 
(MOS) spectra were fitted in the 0.35-15 keV (0.5-10 keV) 
spectral range. 

The Chandra data were reduced with CiAO 3.1, start- 
ing from the linearized "level 2" event lists, using Caldb 
V.2.28. The ACIS-S spectra have been rebinned accord- 
ing to the same criteria employed for the EPIC ones. Fits 
were performed in the 0.1-11 keV energy band. 

Out of the 26 objects in Tab. □ 25 (22) were detected 
at a level larger then "ia in the 0.5-2 keV (2-10 keV) 
energy band. 

4. Spectral results 

Previous studies show that X-ray spectra of obscured 
AGN are complex IjTurner et al. 19971 iRisaliti 2002. 
IGuainazzi et al. 2005a|l . This holds for the objects dis- 
cussed in this paper as well. In Fig. ^ and |21 we show 
the sample source spectra, and their residuals when a 



photoelectrically absorbed power-law model is applied. 
This simple model leaves significant residuals in almost 
all objects. As already discussed by, e.g., Guainazzi et 
al. (2005a), at least two continua components in the 0.5- 
10 keV energy band are required. Fluorescent emission 
lines from iron and other elements, often covering a wide 
range of ionization stages, are often present as well, as 
witnessed by local excesses in the residuals of Fig.^andU] 
at -0.6, 0.9, and 6.4 keV. 

Hence, we have defined a set of physically-motivated 
"baseline" models to fit the data. For Compton-thin 
Seyfert 2s the baseline model is expressed by the following 
general formula: 

FtMniE) = M + + e--(^)^« + ^ G.,{E) 

■i 

where cr{E) is the photoelectric cross-section, with solar 
abundances according to Andres & Grevesse (1989); A is 
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2MASSJD 1500266-0725482 



IRAS154B0-0344 




O.B 1 



1 — 








— ^ 


1 ' — 











- 1 


1 










-1 — 1 1 1- 


— 1 — 1 ; 







=l=d- 













0.5 1 



O.B 1 



Fig. 1. Spectra {upper panels) and residuals against a power-law continuum {lower panels) for the XMM-Newton/EPIC 
observations of our sample. 
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Energy (keV) 

Fig. 2. Spectra [upper panels) and residuals against a power-law continuum (lower panels) for the Chandra/ ACIS 
observations of our sample. 



a normalization factor; fg ("scattering fraction") is con- 
strained: < /s < 1; M is the thermal emission from a 
collisionally ionized plasma with temperature T (we used 
the mekal implementation in XsPEC: iMewe et al. 1985|l : 
Gi{E) are Gaussian profiles, corresponding to fluorescence 
emission lines of high-Z elements from O to Fe, the last 
being by far the most important and common with its 
6.4 keV neutral fluorescent line. For "Compton-thick" ob- 
jects, the baseline model is slightly different: 

FtMck (E) = M + BE-^ + RiE,T) + J2G^iE) 

i 

where i? is a normalization constant, and R{E) is the 
"bare" Compton-reflection of an otherwise invisible nu- 
clear continuum (we used the implementation pexrav in 
XsPEC, [Magdziarz fc Zdziarski 1995| ). In this component, 
we have assumed again solar abundances, and a face-on 
slab. No high-energy cut-off for the intrinsic power-law 
continuum emission has been imposed, as it normally lays 
beyond the sensitive bandpass of the XMM-Newton/EPIC 
or Chandra/ KGIS cameras IjPerola et al. 2002)1 . In princi- 
ple, the Compton reflection component should have been 
added also to the baseline model for Compton-thin ob- 
jects. However, its contribution in the pn energy band- 



pass is expected in these cases to be negligible, and was 
therefore not included for simplicity. 

Observed fluxes in the 0.5-2.0 keV and 2-10 keV en- 
ergy bands, and absorption-corrected luminosity for the 
thermal and the AGN power-law components are shown 
in Tab. 13 

The spectra of the large majority of the sources dis- 
cussed in this paper are comparatively weak, with an over- 
all number of independent spectral channels often lower 
than 50. In most cases, not all the components of the 
baseline model are required. Moreover, in many objects 
it is not possible to discriminate whether the object is 
Compton-thick or -thin from the pure statistical point of 
view. In order to proper classify these objects, we need 
to complement the statistical results of the fits with the 
evaluation of two physical observables with a high diagnos- 
tic power in X-ray obscured AGN: the Equivalent Width 
{EW) of the Fe fluorescent line, and the ratio T be- 
tween the 2-10 keV flux and the Balmer-Decrement cor- 
rected 0[iii] flux. 

The binning requirement imposed by the application 
of the test in the global flts may hamper the detectabil- 
ity of even large EW emission lines in low-statistics spec- 
tra with almost no continuum besides the line. Before ap- 
plying the baseline models to the broadband X-ray spec- 
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Table 2. X-ray observed fluxes, T values, and absorption-corrected luminosities in the 0.5-10 keV energy range 



Source Observed fluxes Luminosities 





0.5-2 keV 


2-10 keV 


5-10 keV'' 


T 


Thermal'' 


AGN'* 


2MASSJ01500266-0725482 


0.30 


0.82 


3.7 


1.00 


9±n K 


1 20±!!'iii? 


ES02nQ-Gni 2 


0.011 


0.06 


4.2 


1.2" 


1.0 ± 0.3 


13±!!,V^ 

u.ionzQ 07 


TGSfiSQ 


0.10 


0.08 


0.39 


0.03 


2 4±?-? 




IC4995 


0.06 


0.29 


2.5 


0.18 


4 ± 2 




IRAS 15480-0344 


0.09 


0.37 


2.8 


0.07 


27±??^ 




Mknl 


0.10 


0.13 


1.0 


0.28 




1 29±?'92 
•^^-'-1 .28 


Mkn266 


0.03 


<0.007 


<0.03 


<0.03 


3 6±^A 




Mkn612 


0.05 


0.36 


3.2 


0.20 


2 7±nf 
^- ' -"-0.5 


7 3±}'? 
' •^-'-1.1 


NGC17 


0.07 


0.33 


2.0 


0.01 


2.0 ± 0.5 


2.16 ± 0.06 


NGC591 


0.09 


0.20 


1.6 


0.11 






NGC1410 


0.04 


0.04 


0.2 


0.01 






NGC3982 


0.018 


<0.05 


<0.2 


< 0.14 


057 + 0014 




NGC4156 


0.13 


0.18 


0.8 






0.414 ± 0.015 


NGG4968 


0.04 


0.15 


1.2 


0.04 






NGC5283 


0.04 


1.46 


10.3 


2.71 




1 7±S'I 

^- ' -"-0.6 


NGG5427 


<0.002 


<0.11 


<0.11 








NGC5953 


0.06 


<0.006 


<0.02 


< 0.01 






NGC7212 


0.09 


0.69 


5.6 


0.22 


8±^ 




NGC7465 


0.24 


4.15 


29.0 






i.o±g:i 


UGC1214 


0.32 


0.28 


1.7 


0.07 






UGC2456 


0.14 


0.36 


2.8 


0.07 






UGC2608 


0.13 


0.14 


0.8 


0.03 


30±?i 




UGC4229 


0.05 


0.22 


1.7 


0.03 






UGC6527 


0.05 


0.43 


4.0 


0.08 


3.2±i;J 


15±^* 


UGC987 


0.48 


1.21 


6.1 


7.24 




1.09 ±0.05 


UM625 


0.03 


0.02 


0.11 


0.004 


2.3±?:i; 




"in units of 10~^^ erg cm~^ s~ 


-1 













''in units of 10 erg cm ^ s ^ 
"in units of 10"'° erg s~^ 
''in units of lO'*^ erg s"^ 

"not corrected for Balmer decrement (measurement not available in the literature) 



tra, we have therefore fit the unbinned spectra in the 
5.25-7.25 keV energy range (where the Ka fluorescent 
iron line complex is expected to lay for all objects in 
our sample), using the Cash statistics IjCash 19761 "lo- 
cal fits" hereinafter). Power-laws were used to indepen- 
dently model the continuum and the background, as the 
Cash method cannot be applied on spectra, which do not 
follow the Poisson statistics. Their parameters were then 
allowed to vary in the local fits within their confidence 
intervals. The properties of the iron Kq, lines derived from 
the local fits are summarized in Tab. |31 Fe Kq emission 
lines are detected in 6 (14) objects at a confidence level 
larger than 2cr (Icr). They represent the first detection 
of such a feature in these objects. In all cases, the cen- 
troid energies Ei are consistent with emission from neu- 
tral or mildly ionized iron. However, in UGC 1214 an ad- 
ditional line with Ei ~6.66 keV suggests a contribution 
from iron species as ionized as Fe xxv. In Fig. O the 
5.25-7.25 keV spectra for the objects where a line is sig- 
nificantly detected are shown. If the Fe EW is larger 
then and inconsistent with 600 eV at the 90% confidence 
level (the maximum observed EW in "Compton-thin" ob- 
jects; EmETeFHIiMZl IEiintr2002|) models assuming 



Compton-thick obscuration have been preferred, whenever 
statistically equivalent to those requiring a Compton-thin 
absorber. 

Another powerful diagnostic of nuclear X-ray obscura- 
tion is the 2-10 keV to 0[ili] flux ratio T. Previous 
studies (jMaiolino e t al. 19981 IBassani et al. 19991 
IGuainazzi et al. 2005a.i) have shown that objects with 
T < 0.1 are invariably Compton-thick, whereas objects 
with T > 1 are almost exclusively Compton-thin or 
unobscured. For all objects where the iron line EW gives 
a constraint on the nature of the nuclear obscuration, 
this is consistent with that inferred from the T value. 
This gives us confidence in using T as a diagnostic of the 
nature of the X-ray obscuration. The values of T for the 
objects of our sample are listed in Tab. HJ 

A summary of the best-fit parameters and results is 
shown in Tab.^ where the column density for "Compton- 
thick" objects is lower-bounded by a^^ ~ 1.6 x lO^'' cm^^. 
The EW, El and intensity // of the Fe Ka fluorescent 
lines derived from the global fits are shown in Tab. 13 
whereas the same properties for other lines are listed in 
Tab. El for the two objects (Mknl and UGC2456) where 
they are statistically requested. We stress that the mean- 



Guainazzi et al.: X-ray obscured AGN in the local universe 

Table 3. Fe K^, properties 



7 





Local fits 




Global fits 




Source 


El 


EW 


El 


EW 






(kcV) 


(keV) 


(keV) 


(kcV) 




2MASSJ01500266-0725482 


6.4''" 


< 200 


6.56±^:1^ 


350^290 


0.23 


ESO269-G012 


6.4+ 


< 700 


6.4+ 


< 3000 


< 0.3 


IC3639 


6.39 ± 0.07 


1500 ± 1100 


6.4+ 


< 525000 


<174 


IC4995 


6.42±S-8? 


1700 ± 700 


6.40±8:8i 


4300±lfgg 


1.1 


IRAS15480-0344 


6.41±S-81 


17000±f8g8o 

1 ■-"-'■-'—1— 


6.4+ 


< 2400 


< 1.5 


Mknl 


6.4+ 


< 800 


6.4+ 


< 2000 


< 0.5 


Mkn266 


6 38±!!-!iI 


6000±2J!S 


6.4+ 


< 100000 


< 50 


Mkn612 


6.46±^Si 


280±?^S 


6.40 ±0.18 


200±?^!] 


0.5 


NGC17 


6.4+ 


< 1900 


6.4+ 


< 322 


< 0.7 


NGC591 


6.41±S-8? 


2200±Inn 


6.47 ± 0.03 


5800±^|8g 


0.9 


NGC1410 


6.420 ± 0.010 


800±I!IS 


6.4+ 


< 11000 


< 0.5 


NGC3982 


6.4+ 


8000 + 5000 


6.4+ 


< 41000 


< 3.0 


NGC4156 


6.4+ 


< 300 


6.4+ 


< 500 


< 0.10 


NGC4968 


6 40±S Si 


3000±}^SS 


6.41 ±0.02 


5300±^^gg 


1.8 


NGC5283 


6.4+ 


< 170 


6.4+ 


< 220 


< 0.7 


NGC5953 












NGC7212 


6.41 ± 0.03 


900±inn 


6.4i±8:8i 


1100 ± 200 


0.9 


NGC7465 


6.4+ 


< 300 


6.4+ 


< 400 


< 3 


UGC1214 


6.41 ±0.05 


1300 ± 500 


6.41 ±0.05 


2700 ± 1000 


0.6 




6.66 ± 0.05 


7oo±!iS!S! 


6.92 ±0.09 


15000±iJ]|]J] 


0.4 


UGC2456 


6.41±°°* 


iooo±tE!8 


6.42±8;!]t 


320 ± 90 


0.6 


UGC2608 


6.4+ 


< 1300 


6.4+ 


< 8000 


< 1.2 


UGC4229 


6.4+ 


< 4000 


6.4+ 


< 1600 


< 0.7 


UGC6527 


6.37±E!:8! 


350±i8 


6.44 ± 0.08 


100 ± 50 


0.3 


UGC987 


6.45 ± 0.07 


180±1J8 


6.48±8:§^ 


120 ± 80 


0.18 


UM625 


6.4+ 


< 1100 


6.4+ 


< 6000 


< 0.2 



+ in units of 10 ® ph cm 
+ fixed 



to NGC4968 




8 6.5 
Energy (keV) 



6 6.5 
Energy (keV) 



6 6.5 
Energy (keV) 



Fig. 3. Zoom of the spectra in the 5.25-7.25 keV (observer's frame) for the sources, where the iron hne is detected 
at a confidence level >2a. The spectra are Unearly binned with a constant SE = 50 eV. Each channel corresponds 
therefore approximately to one- third of a resolution element. 
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Table 4. Global fit continuum spectral results. "G^" refers to the Table number in this paper, where properties of 
the emission lines detected in the global fits are presented. The usage of two rows for a given objects means that two 
thermal components are required by the fit. 



Source Nh 




r 




(10^2 cm-2) 


(%) 


r, 1 U.i3 


(kcV) 



ESO269-G012 
IC3639 
IC3639 
IC4995 

IRAS 15480-0344 
Mknl 
Mkn266 
Mkn612 

NGC1410 

NGC17 

NGC591 

NGC3982 
NGC4156 
NGC4968 
NGC5283 
NGC5953 
NGC7212 

NGC7465 

UGC1214 

UGC2456 
UGC2608 

UGC4229 
UGC6527 
UGC987 
UM625 



>160 

> 160 

> 160 

> 110 

> 160 
65±i 

> 160 
40±^g 

> 160 

> 160 
< 0.02 

> 160 

15 ±3 

> 160 

46±?J 
2.0±[]:i 

> 160 

90±ig 

> 160 

> 160 
76±?^ 

0.068 ± 0.013 

> 160 



o.59±S:?i 



0.9 ±0.5 



i.6±?:l 

9±3 



o.7±!]:| 



3.9±?:! 



2.1 ± 0.9 

2t 

2.4i±j;:i? 



2.40±g:l« 



4.3±8:i 
i.85±g:?g 

2.4 ±0.9 



1.96±8:gi 

2.4±E;i 

2.3 ±0.5 

3.9±y 

i.5±!!:;? 



2.7 ±0.4 

1.97±g:Sg 
2.3±Ji 

2t 

2.3 ±0.3 
1.70 ±0.04 
4.1 ± 0.4 



0.76±g:i 
0.37±g:gg 

i.o±g:l 
o.o96±g:Elil 

r, 17 I 0.07 
U.i '±0.09 



0.78±[]:ii 
0.089±HSI 

o.62±g:E!| 
o.6i±8:i 

0.63±H? 
0.09 ±0.04 
0.69 ±0.05 

o.59±g:ii 



o.i6±!^:Ji 

0.72±H6 



0.16 ±0.03 

o.4±»:|t 



o.6i±8i| 
i.5±g;i 

0.59±g:§? 



o.23±!;:JJ! 



3 

3,5 



0.8/1 
3.3/4 

15.4/11 
4.1/7 

36.0/32 
0.5/1 

32.6/23 

7.2/6 
46.5/40 
29.6/29 

0.7/1 
207.3/179 
23.7/23 
12.5/13 

0.6/1 
97.0/65 

18.0/24 

85.6/72 

60.7/48 
8.4/10 

9.4/7 
15.3/22 
271.8/238 
8.1/6 



^ fixed 



= o.i9±S]:gi 



ing of the EWs corresponding to the "local" and "global" 
fits in Tab.^is not strictly the same. In the "local" fits, the 
EW is calculated against the background-subtracted con- 
tinuum underlying the line profile; in the "global" fits, it 
represent the EW against the proper continuum to which 
it refers, i.e. the absorption corrected primary power-law 
for Compton-thin, and the Compton-reflcction continuum 
for Compton-thick objects. 

4.1. Notes on individual objects 

IC3639, NGC UW, and UM 625: in these objects T < 0.1, 
suggesting Compton-thick obscuration. The EW of the Kq 
iron line is basically unconstrained. However, the spec- 
trum above 2 keV is typical of unobscured AGN, show- 
ing an observed spectral index F ~ 2. The hard X-ray 
luminosity corresponding to the observed fluxes in these 
sources is ~5-7xl0^'' erg s^^. Longer observations would 



be required to ultimately elucidate the nature of the high- 
energy emission in these objects. They can be reconciled 
with Compton-thick obscuration if scattering is domi- 
nated by a hot plasma, or substantially contaminated by 
serendipitous sources in the host galaxies, such as Ultra- 
Luminous X-Ray binaries or bright Supernova Remnants. 
NGC 7465: the application of the baseline model yields a 
very flat intrinsic spectral index: F 0.7, again unusual in 
AGN. A good fit is obtained if one assumes that a "stan- 
dard" (F = 2) AGN emission is covered by a partial cover- 
ing absorber, with column densities Nh.i ~ 5 x 10^^ cm~^, 
and Nh,2 — 2 x 10^^ cm~^, with covering fractions of 
15 ±}!^ % and 85 ±}g %, respectively. 

5. Discussion 

In this Section we combine the objects for which new 
C/ianrfra/XMM-Newton observations have been presented 
in the previous Section, with measurements extracted 
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Table 5. Gi not-iron emission lines detected in the sam- 
ple. The EW are expressed again the scattering continuum 
only. 



Mknl 



UGC2456 



Fe-L, 0.7-0.8 keV 

El (keV) 0.88 ± 0.02 

EW (keV) 160±|!] 

Ii" 2.9 

Ne IX, 0.915 keV 

El (keV) 

EW (keV) 

/i" 

Ne X, 1.02 keV 
El (keV) 
EW (keV) 
/i" 

Fe-L, 1.03-1.15 keV 
El (keV) 1.02 ± 0.02 
EW (keV) 130±i!] 
0.6 

Mg XI, 1.34 keV 
El (keV) 
EW (keV) 

S XIV, 2.00 keV 
El (keV) 
EW (keV) 



0.805 ±0.016 
140 ± 40 
1.1 

0.92 ± 0.08 

230±ig 
1.4 

1.026±HJ9 
180 ± 40 
0.9 

1.15 ±0.02 
100±tg 
0.6 



1.32 ±0.03 
160±TO 
0.46 



1.90 ±0.06 
120±ig 
0.18 



"in units of 10 counts erg s 
'fixed 




-14 -13 -12 -11 -10 

log (2-10 keV flux, cgs) 

Fig. 4. 2-10 keV flux for the objects or our sample 



from the literature on all the other objects of our sam- 
ple. The sample covers a wide range of 2-10 keV fluxes, 
from 3 X 10~" to 6 x 10~^^ erg cm"^ s~^ (Fig- El- 




1000 

K„ Fe line EW (eV) 



2000 



Fig. 5. Normalized distribution functions for the Fe 
line EW in "dusty" and "not dusty" Seyfert 2s. The 
solid bars indicate typical EW values for Compton-thin 
(Nh < 10^3 
Nh < 10^4 cm 
objects (cf. Fig.inj) 



cm 

-2 



Extreme Compton-thin (10^'^ 
), and Compton-thick {Nh > 10^^ cm" 



5.1. Correlation between X-ray absorption and nuclear 
dust content 

In Fig. we show normalized robust distribution func- 
tions'' for the EW of the Kq, iron line - used here as an 
indicator of X-ray obscuration - for "dusty" and "non- 
dusty" Seyfert 2s. The probability that these distributions 
are drawn from the same parent population is 2 x 10~^. 
The distribution weighted mean is 210 ± 80 eV, and 
1000 ± 300 eV for "dusty" and "non dusty" galaxies, re- 
spectively. The latter value is well above the threshold 
separating Compton-thin from -thick objects. This points 
to an average systematic difference in the X-ray absorbing 
column density depending on the nuclear dust content in 
the host galaxy. 

This difference is mainly due to objects with a column 
density 10^^ < Nh < lO^'^ cm"^ ("Extreme Compton- 
ThiN", ECTN, objects hereafter). The ratio R between 
ECTN with and without nuclear dust is Rectn = 1.8 ± 
0.3. For "Compton-ThicK" (CTK) objects the same ratio 
is a factor of three lower {2a level): Rctk = 0.5 ± 0.2. 
The ratio calculated on the whole sample is 0.77. 

M98 define a further group of galaxies with interest- 
ing dust absorption, as those highly inclined spirals with 
extensive dust lanes on one side of their major axis and 
hardly on the other side {"D- [direction] " in M98). 7 ob- 
jects in our sample fall in this class. Although they are 



^ They are histograms of the distribution functions, where 
each measurement is represented by a Gaussian function, 
whose mean is the measurement value and whose deviation 
is the measurement's statistical uncertainty 
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not "nuclear dusty" according to our classification, it is 
interesting to observe that the correlation between X- 
ray obscuration and the presence of dust becomes even 
more extreme if we include them in the "dusty" sample: 
Rectn = 4.5 ± 0.9, RcTK = 1-0 ± 0.4. The ratio calcu- 
lated on the whole sample is 1.56. 

The existence of a correlation between X-ray obscura- 
tion and dust lanes in the HST images, points to a relation 
between the host galaxy gas/dust, and the matter respon- 
sible for X-ray obscuration in the innermost <5 100 pc. As 
already pointed out by GOl, the correlation is probably 
not direct, i.e. the dust lanes seen in high-resolution HST 
images are not directly responsible for the X-ray obscu- 
ration, as in "DI" objects these lanes do not appear to 
cross the line-of-sight to the active nucleus. The connec- 
tion must therefore be of an "ambient" or "evolutionary" 
nature, with dusty nuclear environments statistically fa- 
voring the formation of X-ray obscuring patches and fila- 
ments. 

The most straightforward interpretation of these re- 
sults is that Compton-thin obscuration, primarily in the 
column density range lO'^'^ < Nh < 10^^ cm~^, is due to 
gas associated to the dust lanes. In other words: we mea- 
sure a Compton-thick X-ray obscuration when our line- 
of-sight intercepts pc-scale compact clouds. We will still 
collectively dub them "torus" in the following, although 
their detailed geometry is largely unknown. Whenever this 
does not occur, high-energy photons have still a large like- 
lihood to cross gas associated with the dust lane clouds, 
encompassing the nuclear environment on a larger scale 
and with a larger covering fraction. 

Alternatively, dust lanes and Compton-thick 
"torii" could be mutually exclusive in an AGN. In 
this case, one might expect that the distribution 
of the Mid InfraRed (MIR) luminosity differs be- 
tween "dusty" and "non dusty" Seyferts. The MIR 
Spectral Energy Distribution of obscured AGN is 
in facts dominated by "torus" reprocessing of the 
high-energy continuum ( Alonso-Hcrrcro ct al. 2001', 
[Efstathiou fc Siebenmorgen 2005, ). In Fig. |B1 we show 
such distributions. The MIR flux density was calculated 
from a simple linear combination of 12/i and 25/i IRAS 
measurements: 



S 



MIR 



25 



The distributions are statistically indistinguishable. A 
pc-scale Compton-thick "torus" is therefore most likely 
present in all radio-quiet AGN, in agreement with 
the ubiquity of unresolved iron fluorescent Kq, emis- 
sion lines in their X-ray spectra ( |Page et al. 2003| 
IJimenez-Baiion et al. 2005j) . 

The viability of dust lanes as tracer of Compton- 
thin X-ray obscuring gas depends on whether 
galactic dust outside the nucleus can produce 
enough extinction. Estimates from HST color 
maps dMulchaey et al. 1994[ [Simpson et al. 1997| 
IFerruit et al. 20001 IQuillen et al. 1999aD suggest 
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6 

4 - 
2 

_ 
10 



4 - 
2 



Dust 




No dust 



42 43 44 45 

log(L„|R, erg s"') 



46 



Fig. 6. MIR luminosity distribution histogram for "dusty" 
{upper panel) and "non-dusty" {lower panel) objects in our 
sample. Shaded areas indicate upper limits. 



i0.5~1.5 mag. This corresponds to Nh 



iO.5- 



1.0x10 cm for gas-to-dust ratios typically 
observed in luminous and obscured nearby AGN 
(IGranato et al. 19971 IMaiol ino et al. 2001|l ^. The fact 
that a large fraction of "dusty" Seyfert 2 in M98 have a 
rather disturbed morphology is suggestive of dumpiness 
of the absorbing medium, which may substantially 
increase the line-of-sight column density in individual 
objects. 

Maiolino et al. (1995) suggested the presence of a 
"100 pc torus" in intermediate {i.e. 1.8-1.9) type Seyferts 
to explain the fact they they are preferentially found in 
edge-on hosts. Following the results of their HST survey, 
M98 elaborated the "Galactic Dust Model" (GDM) for 
obscuration - as opposed to the standard "Torus Model' 
- where the whole obscuration in type 2 AGN is ex- 
plained by the 100 pc dust lanes imaged by HST. However, 
M98 themselves realized that this model is not applica- 
ble to Compton-thick objects, as dynamical considera- 
tions prevents the host galaxy dust lanes from providing 
such extreme obscuration. The GDM is therefore insuffi- 
cient to explain the large fraction of local obscured AGN 
l|Risaliti et al. 19991 this paper). This motivated Matt 
(2000) to propose an extension of the standard unified 
scenario, whereby the compact, dusty "torus" intercepts 
our line-of-sight to Compton-thick AGN only, whereas 
Compton-thin obscuration could be due to gas and dust 
located at larger scales, possible related to the host galaxy 
rather than to the nuclear engine. [Guainazzi et al. 2005al 



^ Interestingly, the average surface density of molecular gas 
in the innermost 500 pc of the Milky Way corresponds to Nh — 



4 X 10^^ cm"^ ISanders et al. 1984J 



Guainazzi et al.: X-ray obscured AGN in the local universe 



11 



and IMatt et al. 2003BI have shown that some Compton- 
thick objects would appear to us as Compton-thin, if we 
could observe them through a different Une-of-sight. This 
supports the presence of a multiplicity of obscuring sys- 
tems in the innermost hundreds parsecs. 

It must be emphasized that the existence of ex- 
tended obscuring structures on ~100 pc scale as such 
does not rule out the existence of an obscuring torus. 
Compact toroidal geometries l|Pier fc Krolik 1993)l 
are unable to explain the broadness of the spectrum 
in the 1-20/Ltm region, and the nature of the 9.7fim 
silicate emission. However, several different combina- 
tions of geometry and torus physical properties are 
consistent with the observed IR SED of highly ob- 
scured AGN IjClalhano et al. 20T)8| . once allowance is 
made for "extended torii" l|Granato fc Danese 1994(1 . 
where the ratio between the outer and the inner 
radius is ~10-100. Smooth HGranato et al. 1997)l 
or clumpy IjNenkova et al. 2002|l density distribu- 
tions, anisotropic flares and tapered disk geometries 
l|Efstathiou fc Rowan-Robinson 1995|l are all equally 
successful in reproducing the data. Moreover, dust clouds 
might be associated as well with the extended ionizing 
cones. Mid-infrared emission is slightly extended in many 
Seyfert 2 IjMaiolino fc Rieke 19951 IGalhano et al. 2003|l . 
Dust lanes are sometimes associated with 0[iii] ionization 
cones as well dQuillen et al. 199*931 ). Models including 
an optically thick torus and an optically thin "cone" 
of dust clouds (Efstathiou & Rowan- Robinson 1995) 
can as well reproduce the observed IR SEDs 

onso-Herrero et al. 2003|) . Finally, radiation pressure 
driven starburst disks could contribute as well to the X-ray 
obscuration on scales ^ 100 pc ( [Thompson et al. 2005| ). 
All these possibilities can hardly be discriminated even 
with the highest resolution imaging detectors currently 
available, but could be disentangled spectroscopically 
l|Lutz et al. 2004| and references therein). 



5.2. Overall sample properties 



The parent sample, from which our sample discussed is 
extracted, was selected on the basis of observational prop- 
erties, which should not introduce any X-ray selection ef- 
fects. Our XMM-Newton program, from which most of 
the XMM-Newton observations discussed in this paper are 
drawn, carefully avoided any possible X-ray related bias. 
However, the inclusion of archival Chandra and XMM- 
Newton observations, belonging to different projects, as 
well as of the GDI objects, undoubtedly may introduce a 
potential bias toward X-ray brighter, and therefore poten- 
tially less obscured, objects. With this caveat in mind, we 
review in the following Section the main X-ray spectral 
properties of our sample. 



1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 


1 1 1 1 1 1 1 


- 





































39.5 40.0 40.5 41.0 41.5 42.0 

log (Lo.3-2kev. erg/s) 

Fig. 7. Distribution function of the 0.3-2 keV intrinsic lu- 
minosity for the sub-sample of our sample, observed by 
either Chandra or XMM-Newton. The filled area corre- 
spond to objects, whose hard X-ray emission is covered 
by a column density Nh < 1.5 x 10^^ cm~^, implying 
contamination by the transmitted nuclear emission in the 
soft X-ray band. 

5.2.1. Soft excess 

The hard X-ray continuum of type 2 Seyferts is dominated 
by a power-law continuum, partly suppressed by photo- 
electric absorption IjAwaki et al. 19911 ITurner et al. 19971 
IRisaliti 2002|l . In Compton-thick Seyferts, the transmit- 
ted continuum is totally suppressed below 10 keV, and 
the hard X-ray continuum is accounted for by reflec- 
tion of the primary nuclear emission off the inner walls 
of optically thick matter surrounding the nuclear en- 
vironment (^Matt 2000 ). The spectra of our sample - 
including those objects whose X-ray spectra are pre- 
sented for the first time in this paper - are consis- 
tent with this scenario. Highly-obscured AGN exhibit al- 
most invariably large soft excess emission below 2 keV 
IjGuainazzi et al. 20 05a'). In our sample, the total inte- 
grated soft X-ray luminosity of this component range be- 
tween 10**° erg s~^ and 4 x lO'*^ erg s~^ (Fig. CJ. In prin- 
ciple, this excess emission could be associated with un- 
resolved sources or diffuse emission in the host galaxy, 
outshone by the formidable output of the AGN in unob- 
scured objects. Although still tenable for some individ- 
ual objects, this explanation cannot apply to the sam- 
ple as a whole. Chandra measurements of the luminos- 
ity function in nearby spirals suggest that at most a few 
sources with Lx ^ 10^° erg s""'^ are present in each galaxy 
HBauer et al. 20011 see lWolter fc Trinchieri 2004l for an in- 
triguing exception: the Cartwheel galaxy). Likewise, there 
is little - if any - diffuse emission associated with spiral 
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NGC5Z83 - Log(LX)~40-Z NGC5953 - Log(LX)~40.2 Mkn256 - Log(LX)~40.8 

HST (greyscale) — Chandra 0.1-2 keV (contours) H5T (greyscale) — Chandra 0.1—2 keV (contours) H5T (greyscale) — Chandra 0.1-2 keV (contours) 




3 29 55 81 100 4.7 7.5 11 21 100 3 29 55 81 100 

Fig. 8. HST FW606 filter (greyscale) and Chandra 0.1-2 keV images (contours) for: NGC 5283 (left), NGC 5953 
(center), Mkn 266 (right). Chandra images were smoothed with a 0.75" Gaussian kernel. Contours levels are linearly 
spaced from 1 to the maximum of measured pixel counts. 



galaxies (at the level of a few lO'^^ erg s ^ in M 31; 
LPictsch ct al. 2005, Trinchicri 2004). We have searched 
for diffuse soft X-ray emission in the 7 objects of our 
sample observed with Chandra^. In 4 objects the radial 
profile is consistent with the instrumental PSF. In the 
two galaxies with the lowest 0.3-2 keV luminosity (1.5- 
1.7x10^° erg s^^) we indeed detect a large-scale low- 
surface brightness component (left ad central panel of 
Fig. inj, that in NGC 5953 extends across a significant 
fraction of the optical surface of the host galaxy. Finally, 
a peculiar object is Mkn 266 (right panel of Fig.|Sl), which 
belongs to a close interacting galaxy pair, embedded in a 
luminous (7x10'*^" erg s""'^) X-ray halo. A detailed analysis 
of this intriguing system is deferred to a future paper. 

In a few bright, nearby Seyfert 2s, for which 
high resolution spectroscopic measurements have 
been possible (|Sako et al. 20001 ISambruna et al. 20011 
IKinkhabwala et al. 2002|l . soft X-rays are dominated by 
emission lines from highly ionized species due mainly 
to recombination and resonant scattering following 
photoionization and photoexcitation. The density of the 
plasma is generally not large enough to substantially 
contribute with a scattering continuum. The striking 
coincidence between diffuse soft X-ray and 0[iii] emis- 
sion is consistent with this scenario (Bianchi et al., in 
preparation). However, in some low-luminosity Seyfert 2s 
(Lx < 10^" erg s~^), high-resolution imaging observa- 
tions with Chandra have shown that a contribution to 
the soft X-ray by plasma coUisionally ionized in episodes 
of intense nuclear star formation is dominant in the soft 
X-ray regime Ij Jimenez-Bailon et al. 2003)l . It is inter- 
esting to observe that this soft component is generally 



unobscured. It must be therefore produced in or by gas, 
extending well beyond the material absorbing the nuclear 
continuum. This unfortunately provides little constraint 
on its origin, as both nuclear starbursts and AGN-driven 
photoionized extended structures on scales of the order 
of several hundreds parsecs are common in nearby AGN 
((Wilson et al. 1992J ISako et al. 20(101 |Young et al. 200 1[ 
pBianchi ct al. 2003^. 

This ubiquitous soft excess may have an important im- 
pact on the classification of Compton-thin/-thick objects 
at cosmological redshift, if one uses standard X-ray survey 
colors. This point is discussed in Appendix A. 

5.2.2. Ka iron lines 

It has been early recognized that fluorescent Kq, iron lines 
are a very common feature in the spectrum of highly ob- 
scured AGN IjAwaki et al. 199l)l . In Fig. El we present the 
EW of this feature (calculated against the total contin- 
uum) versus the X-ray column densities for all the ob- 
jects of our sample, which do not have a lower limit on 
their column density, as in most cases this lower limit was 
inferred on the basis of the high iron line EW - as dis- 
cussed in Sect. 3. The two quantities are well correlated 
above a column density NH,tr ^ 10^'^ cm^^. The correla- 
tion flattens to an almost constant EWo for Nh < NH,tr- 
The measured EWs are significantly larger than expected 
in transmission from a uniform shell of material encom- 
passing the continuum source ( ILeahy fc Creighton 1993) ). 
On the other hand, obscuration of continuum photons at 
~6.4 keV starts becoming significant at column densities 
— Nn^tr- A censored data fit with a function: 



The typical Point Spread Function (PSF) of the XMM- 
Newton observation correspond to >25 kpc radius 



EW(Nh) = EWoe'^'^-^'''^" 
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Fig. 9. X-ray obscuring column density against the iron 
Kq fluorescent emission line EW (details in text). The 
dashed line shows the predicted EW from a uniform 
shell of material encompassing the continuum source 
jLeahy &: Creighton 1993) ). The dot-dashed line is the 
best-fit in a scenario, where the line is generated by re- 
flection off the inner wall of optically thick matter seen 
along an unobscured line-of-sight, while the underlying 
continuum is obscured by matter with a column density 



where (T6.4 is the photoelectric cross-section at 6.4 keV, 
yields an excellent fit: x^/^ — 17.4/20. The func- 
tion levels off at EWq = 113 ± 13 keV. This sup- 
ports an origin of the iron line from the inner wall of 
the "torus" (Ghisel hni et al. 19941 iKrolik et al. 1994|l . It 
is interesting to observe that EWq is in good agree- 
ment with the EW of the iron Kq, line observed in 
low-luminosity {Lx ^ 10^^ erg s~^), X-ray unob- 
scured quasars Ij Jimenez-Bailon et al. 2005)l . This sug- 
gests a common origin of this feature in X-ray obscured 
and unobscured AGN (jMatt 2000,) . 



5.2.3. Nh distribution 

Compton-thick objects largely dominate the sample below 
a 2-10 keV flux of 5 x 10"^^ erg cm~^ s~^. There is, how- 
ever, no explicit dependency of Nh on the intrinsic AGN 
luminosity. If we use the 0[iii] flux as a proxy for the latter 
quantity (it is impossible to have an estimate of the intrin- 
sic AGN luminosity from X-ray refiection-dominated spec- 
tra) , no correlation is observed with the measured column 
density (Fig.CH see also' Risaliti et al. 1999|l . in agreement 
with the 0th order unified models. 
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Fig. 11. 0[iii] luminosity versus column density for the 
objects of our sample. 



The total fraction of Compton-thick objects in our 
Seyfert 2 sample is 46 ± 10%. This is likely to represent 
a lower limit to the fraction on Compton-thick objects in 
the parent sample, due to our possible bias toward X-ray 
brighter, and therefore potentially less obscured, objects. 

6. Conclusions 

In this paper, we have presented new low-resolution X- 
ray spectroscopic measurements for 26 nearby (z < 0.035) 
Seyfert 2 galaxies. We join these results with information 
extracted from the literature to build up a sample of 49 
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objects. The main goal of this paper is the study of the 
correlation between the X-ray column density and the dust 
content on <, 100 pc scale as imaged by the HST/WFC2. 
Our conclusions can be summarized as follows. 

— the presence of "dust lanes" is correlated with X-ray 
obscuration primarily in the column density range 
10^3 < < 10^'* cm^^. We interpret this results 
as due to the larger covering fraction of the gas 
associated with the dust lanes, with respect to the 
compact, pc-scale, Compton-thick "torus". Objects, 
whose line-of-sight to the nucleus does not intercept 
the "torus", still have a non- negligible chance to be 
absorbed by gas on a larger scale 

— the fraction of Compton-thick objects in our sample is 
46±10% and is strongly (observed-) flux dependent, as 
expected. No correlation exists between the amount of 
X-ray obscuration and isotropic indicators of intrinsic 
AGN power in the luminosity range covered by our 
sample {Lx ^ 10**^ erg s^^) 

— the spectrum of highly obscured AGN is complex, with 
a prominent soft excess above the absorbed/reflected 
nuclear continuum. This excess largely dominates the 
overall X-ray observed energy output. This soft compo- 
nent is generally unobscured. This suggests an origin 
in gas structures, which extend beyond the material re- 
sponsible for the obscuration of the nuclear continuum. 
The presence of this soft excess may confuse the esti- 
mate of the column density in highly-obscured, z <, 1 
objects, when standard X-ray colors are used, as shown 
in Appendix A. 
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Appendix A: Brief note on the usage of X-ray colors 
to recognize Compton-thick objects 

Regardless of the detailed physical origin of the soft ex- 
cess, it is interesting to estimate how such a component 
- which can account for a significant fraction of the ob- 
served flux in the X-ray band - may affect the identifi- 
cation of highly obscured AGN at z ^ 1 (at larger red- 
shifts this component is shifted outside the energy bands 



where X-ray surveys are conducted). In Fig. 1121 we show 
color plots for the sample of Compton-thick AGN discov- 
ered by the analysis described in our sample (Sect. 4) to- 
gether with the optically-defined sample of Compton-thick 
Seyfert 2 galaxies discussed by Guainazzi et al. (2005a; the 
Circinus Galaxy has been excluded from the sample, as its 
soft X-ray emission is substantially suppressed by a large 
Milky Way's intervening column density). For each source, 
we have measured the count rates in 4 different energy 
bands: 0.2-0.5 keV (Ultra-Soft; US), 0.5-2.0 keV (Soft; 
S), 2.0-4.5 keV (Medium, M), and 4.5-10.0 keV (Hard, 
M). We have calculated the following Hardness Ratios 
(HRs: |Hasinger et al. 200il ): HRl = {S -US)/{S + US), 
HR2 = {M- S)/{M + S), Hm = {H ~ M)/{H + M). 
Fig.ll2lshows the soft color plane {HR2 versus HRl) and 
the hard color plane {HR3 versus HR2) for Compton- 
thick "local" galaxies, together with the loci occupied by 
"type 1" AGN, "type 2" AGN and unclassified highly ob- 
scured (logiNn) > lO^i-^) sources in the XMM-Newton 
observation of the Lockman Hole IjMainieri et al. 2008|l . 
Compton-thick Seyfert 2 galaxies occupy a specific area 
of these planes, thanks to their combination of large HR3 
(flat hard X-ray spectrum) and low HR2 (steep soft X-ray 
excess). Interestingly enough, they span a range in HRl 
which encompasses typical values for both "type 1" and 
"type 2" AGN, suggesting a multiplicity of possible origins 
for the soft excess emission. In the same Figure we show 
which loci our "local" Compton-thick sources would oc- 
cupy, if they were located at z = 1. In the hard color plane, 
Compton-thick Seyfert 2 galaxies would be basically indis- 
tinguishable from other AGN types. A sizable fraction of 
them would actually occupy the same region as broad-line 
AGN. The soft color plane would be required to remove 
the degeneracy, although Compton-thick Seyfert 2s would 
still have a softer HR2 with respect to other less obscured 
AGN, which could lead to an underestimate of the column 
density covering their active nucleus, and consequently of 
the fraction of highly obscured AGN at z ~ 1 . 
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